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Abstract 

In this paper we study a so-called separatrix map introduced by Zaslav- 
skii-Filonenko [ZF68| and studied by Treschev and Piftankin |Tre981[Trer)21 
[plfnHlfpTnT] . We derive a second order expansion of this map for trigono¬ 
metric perturbations. In |CK15|lGK15| . and |KZZ15j . applying the results 
of the present paper, we describe a class of nearly integrable deterministic 
systems with stochastic diffusive behavior. 
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1 Introduction 


The main goal of this paper is to derive a second order expansion of a so-called 
separatrix map for a class of nearly integrable systems. In nearly integrable 
Hamiltonian systems with one and a half degree of freedo m this map was in¬ 
troduced by Zaslavskii and Filonenko in |ZF68] . Shilnikov Sil65 , using a very 


similar geometric idea, studied it in a neighborhood of homoclinic orbits without 
restriction to Hamiltonian structure or closeness to integrability. Treschev and 
Piftankin estimated the error terms in the traditional version of the separatrix 
map and studied the multidimensional situation |Tre98[ ITre02[ IPif06| . 

It is becoming clear that the separatrix map is a powerful tool to analyze 
the dynamics in a neighbourhood of homoclinic orbits to a normally hyperbolic 
invariant manifold and instabilites of Hamiltonian systems (see the survey prnT] 
and the papers by Treschev |Tre04l ITrel2] b For this reason in |KZZ15| . using 
results of this paper, we perform an indepth analysis of the phenomenon of global 
instabilities in nearly integrable Hamiltonian systems. Usually this phenomenon, 
discovered by Arnold |Arn64] . is called Arnold diffusion. 

The purpose of this paper is to have detailed studies of the multidimensional 
separatrix map, proposed in |Tre02] . and to compute higher order expansion of 
this map with smaller remainder terms. 

The main motivation for this work is to study certain stochastic diffusive 
behavior for nearly integrable deterministic systems. Such stochastic behavior 
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was conjectured by Chirikov |Chi79] in the 1970s. 

Recall a basic fact from probability theory, which is a non-homogenenous 
version of Donsker’s theorem (see e.g. |Eth05] ). Let r/ G M, a{ri) > 0 and 6 ( 77 ) be 
two smooth functions of one variable. Fix 5 > 0, t]o and consider a sequence 

Vn+l = Vn + (j{r]n)SUn + 

where UnS are independent random variables taking values 1 or —1 with equal 
probabilities. Then for each s > 0 and ns = as 5 —)■ 0 the distribution of 

ring converges weakly to the distribution of the Ito diffusion process starting at 
rjo with the drift b{r]) and the variance cr^(? 7 ). 

In the framework that we study in this paper, rj is an action variable of some 
Hamiltonian system. To study long time evolution of r]t we consider a separatrix 
map and study a discrete version rjt^, n > 1. Then the variance < 7 ( 77 ) of the 
discrete version is given by an associated Melnikov function and was computed 
by Zaslavski (see (|^ for a more precise claim). However, to analyze rjt as a 
diffusion process we also need to compute the drift function 6 ( 77 ). This is the 
main purpose of our paper. 

In Appendix we apply our analysis to the so-called generalized Arnold ex¬ 
ample and compute analogs of the drift 6 ( 77 ) and the variance < 7 ( 77 ). In |KZZ15j 
we combine these results with results from |CK15j (see also |GK15j l and con¬ 
struct a class of probability measures for a class of trigonometric perturbations of 
the generalized Arnold example. It turns out that the distributions of a certain 
action component under the associated Hamiltonian flow with respect to these 
measures weakly converge to a stochastic diffusion process on the line. 

More precisely, in |KZZ15] for an open class of trigonometric perturbations 
for a time one map of a generalized Arnold example, which is a 4-dimensional 
symplectic map we construct a normally hyperbolic invariant lamination A. 
Leaves of this lamination are diffeomorphic to 2-dimension cylinders A = T x 
M. We hnd a coordinate change such that the restriction oi F = to this 
lamination is a skew-product of the form 

F:Ax{0,1}^, F{x,u:) = {f^{x),(Tu:), 

where u G {0,1}^, <7 is the shift, and fuj{x) is an exact area-preserving map of 
the 2-dimensional cylinder A. Loosely speaking, using the results of the present 
paper, we obtain a second order expansion of fjs in the perturbation parameter. 

In |CK15j and |GK15j we analyze a class of skew products, which includes 
the aforementioned system F{x,oj), and prove weak convergence to a stochastic 
diffusion process on the lin^ 

^Examples of nearly intergrable systems with stochastic diffusive behavior were constructed 
by Marco and Sauzin (see [MSOdl [Sau06| L 
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We emphasize that the previously known results on Arnold diffusion for the 
generalized Arnold example or, more generally, for a priori chaotic systems, or 
a priori unstable systems, or a priori stable systems, establish the existence of 


“special diffusing” orbits (see [RernSl IRK^ IRTM IRT^ ICYM ICYTM \C7A^ 
Klheldl IDdlLSnni IDdlLSObi [DM IDdlLS13i IFClKPlll iGTiM IGdlLObl 101041 
IKain.Sl IKMVM IKT.nSal iKT.nSbl IKST^ iKm [irhSTdl IMSM IMat91al IMatQlbl 
IMat931 [Mat961 [Mat031 iMatOSlIMoe961 IPifOGl ITre041 ITrel2] 1. In |KZZ15j . |CK15j . 
and |GK15j . using the analysis of the present paper, we study deterministic evolu¬ 
tion of a class of probability measures under the deperministic Hamiltonian flow 
and obtain that the distributions of a certain action component weakly converge 
to a stochastic diffusion process. 


1.1 The separatrix map for the Arnold example 


Since the definition and formulas of the separatrix map are rather involved, we 
start by considering a simple case: the generalized Arnold example. Later, in 
Section 1^, we consider more general models. In this section we consider Hamil¬ 
tonians of the form Hq + sHi with 




J2 2 

y + y + (cOSg-l), 
P(/,p,e^e*^e‘0, 


( 1 ) 


where the unperturbed Hq defines the product of the rotor and the pendulum 
and P(/,p, e*^, e*‘^, e®*) is a real valued trigonometric polynomial in ip and t. For 
the classical Arnold example P = (1 — cosq){cos(p -|- cost). 

In this section, we describe and provide “simple” formulas for the separatrix 
map for such models. The rigorous definition is provided in Section and more 
accurate formulas for model ([^ are derived in Appendix 

The Hamiltonian of the pendulum y -|- (cos g — 1) has a saddle at (p, q) = (0,0) 
whose stable and unstable invariant manifolds coincide forming a figure eight 
(see left picture in Figure [^. In the full phase space, the saddle corresponds to 
the normally hyperbolic invariant manifold {p,q) = (0,0), I G [/_,/+], G T. 
The union of the figure eights for every value of / G [/_, /+] form the stable and 
unstable manifolds of the normally hyperbolic cylinder. Consider a neighborhood 
of such set. It can be defined by 



-|- (cos g — 1) 


< c. 


/ G [/_,/+],(p G T 


for some small c > 0. When one adds the perturbation Hi, typically, the invariant 
manifolds of the cylinder split. The goal of the separatrix map is to understand 
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Figure 1: The figure 8. 


the dynamics in the neighborhood of the perturbed invariant manifolds, which 
now intersect transversally. 

Now we describe the separatrix map. A more rigorous definition of the sep- 
aratrix map is done in Section Consider the time one map associated to 
the flow of the Hamiltonian H^. We define two fundamental domains, as shown 
in Figure]^ as domains such that the image of any point in them under does 
not belong to them and some iterate of any point near an unstable manifold in¬ 
tersects one of these domains. For any point in the fundamental domains one 
can keep iterating T^. If some iterate of such point enters again into one of the 
fundamental domains, we call this point the image of the initial point by the 
separatrix map. Note that some points may not have any of the further iterates 
in the fundamental domains, e.g. those belonging to stable invariant manifolds. 
Then, for these points, we say that the separatrix map is not defined. Sometimes, 
a map of this kind is also called an induced map or a return map. 

To give formulas for the separatrix map, we consider the Fourier expansion of 
the Hamiltonian Hi{I, (p, 0, 0, t), 


\k\<N 

where > 0 is the degree of the trigonometric polynomial. 

Let 'ijj : M —> [0,1] be a function such that '0(r) = 0 for any |r| > 1 and 
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Figure 2: The fundamental domain 


'0(r) = 1 for any |r| < 1/2. We define 


|fc|<Ar VP/ 


for some fixed /? > 0 independent of e small enough so that for every I G [/_, /+] 
only one harmonic does not vanish (we do not have overlapping of resonant zones). 

We consider a certain system of coordinates [r], h, r, a) in the fundamental 
domain (it is explained more precisely in Theorem 2.1). The last variable a G 
{ —, +} denotes the two connected components of the fundamental domain. We 
define the function 


Wo{r]*,h*,C) = h*- 






Then, in such coordinates, the separatrix map is dehned for points such that 
kol e and has the following (implicit) form 


T]* =1] 

r=e 

h*=h 

T* =T 

a* = a sgn wq 


+ ed^e^{T], r) + wo, r) + 

+ rjloglK^^wol + 0{e\oge) 

+ edrQ^iT],^,T) + 0{e^/^) 

+ \og\K'"wo\ + 0{e\oge) 
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where Q are certain constants (see Section 2.1 and Appendix [A]), 


i + r]log\K^Wo\)) 


and are the Melnikov potentials, which are defined as 



©'"(h, ^,r) = / (i^i (T^iv, + - Hi{r],(p + r]t, 0, 0, t - r)) dt, 


where are the time parameterization of the pendnlnm separatrices, that is 



Now we consider a more general set np, proposed by Treschev. 


1.2 The set up 

Consider a Hamiltonian system 




( 2 ) 


where / G M" are actions, (p G are angles and {p, q) belong to an open domain 
D C Even if not written explicitly, the Hamiltonian Hi may depend on the 
parameter £. Fix a bonnded open set T> C M”. 

We assnme that for every fixed / G "D, the Hamiltonian HQ{I,p,q) has a 
saddle at {p,q) = (0,0) with two separatrix loops (see Fig. [^. In |Tre02j it is 
assnmed the following hypotheses. 

HI The fnnction H is C^-smooth in all argnments while Hq is real-analytic in 
p, q, and C^-smooth in I. 

We consider the alternative assnmption. 

HI' The fnnction Hq is and H is C®-smooth in all argnments for s > 6 and 
r > 8s -1- 2. 

That is, we admit lower regnlarity on Hq. On Hi we assnme one degree 
more of regnlarity than in |Tre02] . It is needed to have better estimates of the 
separatrix map. 
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H2 For all points G P, the function H{I^,p,q) has a non-degenerate saddle 
point at (p, q) = (p°, g°) smoothly depending on I. For all G V, (p°, g°) 
belongs to a connected component of the set 

{(p,g): i/o(Ap,g)=^o(/°,P°,g°)}. 

Moreover, is the unique critical point of Ho{I^,p,q) in this compo¬ 

nent. 

Remark 1.1. Using Prop.l, \Tre02j . if one assumes that the saddle is at a certain 
point (p, q) = (p°, g°) which depends smoothly on I, then, one can perform a 
symplectic change of coordinates so that the critical point is at (p, q) = (0, 0) for 
all I eT>. After such a coordinate change in HI is replaced by 

We denote the loops of the “eight” given by Hypothesis H2 by 7 ^(/°). These 
loops have the natural orientation generated by the flow of the system (see Fig. 
[^. We can define an orientation in D by the coordinate system {p,q). 

H3 For all 1° G P, the natural orientation of 7 ^(/°) coincides with the orienta¬ 
tion of the domain, i. e. the motion of the separatrices is counterclockwise. 

In |Tren2] . Treschev defines the separatrix map for Hamiltonians satisfying 
these hypotheses and obtains a formula for this map with certain remainder 
terms. The goal of this paper is to refine Treschev formulas in several aspects. 

Note that the formulas for the separatrix map present certain differences in 
what are called a non-resonant and a resonant regime (see below). Treschev 
provides global formulas for the separatrix map. Here we separate the two regimes 
and give more precise formulas. The refinements we do are the following. 

• For the non-resonant regime we compute the separatrix map up to 2nd 
order in e. 

• For the resonant regime we give formulas in slow-fast variables and we 
improve the size of the remainders. 

We obtain formulas for general Hamiltonians, but also pay attention to the par¬ 
ticular case of the generalized Arnold example ([^. The results for such model 
are presented in Appendix [A] 

Acknowledgement The authors thank Dmitry Treschev for helpful discus¬ 
sions and remarks on the preliminary version of the paper. The authors also 
acknowledge many useful discussion with Ke Zhang. The first author is partially 
supported by the Spanish MINECO-FEDER Grant MTM2012-31714 and the 
Catalan Grant 2014SGR504. The second author acknowledges NSF for partial 
support grant DMS-5237860. 






2 The separatrix map: Treschev’s results 

We devote this section to define the separatrix map and state the results obtained 
in |Tren2] . 

We want to define a separatrix map for points whose {p, q')-components are 
“near” the unperturbed separatrices (see the shaded region on Fig. below). 
For the unperturbed system there is a normally hyperbolic invariant cylinder 
A = {p = q = 0} and for each 1° G P we have invariant tori 

A(/0) ;= {/ = = g = 0} = A n {/ = 1°}. 

These tori are (partially) hyperbolic, i.e there are expanding and contracting 
directions dominating the other ones. There exist two asymptotic manifolds, 

f=^(/°) C {{I^,ip,p,q,t) : ip,teT, Ho{I^,p,q) = 0} 
f±(/0) = {/0}xTx7±(/o)xT, 

where 7^(/o) cire the two separatrices in the {p,q) plane. The manifolds F^(/‘^) 
consist of unperturbed solutions that approach A(J°). 

Assume that P is an open connected domain with compact closure P. In 
what follows, we consider the dynamics of the non-perturbed system in a neigh¬ 
bourhood of the unstable and stable manifolds of the cylinder A: 

r = u „5 (f+(/)uf-(/)). 

This neighbourhood contains the most interesting part of the perturbed dynamics. 
It is convenient to pass to the time one map T^: for any point (/, 93 , p, q) 

Te : {I,(p,p,q) —^ (/(l),(p(l),p(l),g(l)), 

where {I(t), ip(t), p(t), q(t)) is the solution of the Hamiltonian system given by 
He with initial conditions (/(O), (p(0),p(0), g(0)) = {I,ip,p,q). The map Tq has 
1-dimensional hyperbolic tori L{I) = 7r(A(J)), where the map tt : {I,ip,p,q,t) —?• 
{I,ip,p,q) is the natural projection. Let F1=*"(J) = 7r(F^(/)) be the invariant 
manifolds (see Fig. [^. We define the separatrix map SAA^ corresponding to 
in a neighbourhood of the set 

E = u„5(E+(/)uE-(a) = Xr)- 

Let 17 be a small neighbourhood of the set Ujg^L(/) and let U be a neigh¬ 
bourhood of S (see Figure where S is pictured in the case n = 0). If U is 
sufficiently small, then U \ 17 consists of two connected components and U^, 
and, thus, C U U. 
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Figure 3: Neighbourhoods about the hgure 8. 

Consider a point z G U+UU~. Let ki = ki{z) be the minimal positive integer 
such that ^ U U~ and let k 2 = k 2 {z) be the minimal positive integer 

such that k 2 > ki and T^^{z) G U+ U U^. The trajectory T^{z) leaves U+ U 
at k = ki, and it returns to U+ U at k = k 2 - A point z is said to be good if 
k 2 < +C )0 and ... ,T^'^~^{z) G U. Setting 

Ue = {z G U U“ : z is good} 

we obtain maps 

SM,i-,k2i-) + k) : U, ^U+UU~, 

SMe(z, k2{z) + k)= 

Here k G {0,1,... } is a parameter, and it is assumed that 
2^fc2(d+i(^)^... G U+ U U-. 

In |PTn7j it is set t+ = k 2 {z)+k. Since neither U"'' nor U~ serves as a fundamental 
domain, there is a considerable freedom for k. We would like to avoid this freedom. 
By analogy with the pendulum case (see Fig. we let A"*" and A“ be hyperplanes 
in A X A X T 3 (J, (p, p, g, t) whose projection onto the (p, g)-component are curves 
going from one of the connected components of the boundary of |iLo| < c to 
another and transversal to the upper and lower separatrix loops. We denote by 
the subdomain of Uc between the curves A^ and T£(A^) (see Fig. |^. We 
choose such k above that G A+ U A“. 

To provide formulas for the separatrix map we need to set up some notation. 
We dehne 

E{I) = Ho{1, 0, 0), iy{I) = diE{I) ; P ^ M” (3) 
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The Hamiltonian Hi{I, (p, 0, 0, t) has the Fourier expansion 


fceZ"+i 


Let "0 : M —> [0,1] be a C°° function such that 0(r) = 0 for any |r| > 1 and 
0(r) = 1 for any |r| < 1/2. We dehne 


fceZ"+i V / 


(4) 


for some constant 0 > 0. 

To obtain quantitative estimates, we follow |Tre02] and use skew norms, de- 
hned as follows. Let K C M”* be a compact set and j G N. Then, for functions 
f e {V X K) we dehne 


11/(1, ^)||f= max 


Ql'+l" 


■ ■ ■ d^'^Ind’-'izi ■ ■ ■ d’-'^Zr, 


where I" = + ■ ■ ■ + /(/• It is assumed that / can take values in M'^, where s is an 

arbitrary positive integer. The norms || ■ ||r^^ are anisotropic, and the variables r 
play a special role in these norms because the additional factor b corresponds to 
the derivatives with respect to r. Obviously, || ■ ||r^^ is the usual C'’-norm. This 
norm is similar to the skew-symmetric norm introduced in |KZ12j . Section 7.2. 

For a function / G C''(P x K) and g G C^{'D x K) we say that 

f = 0<P(g) if ||/||?><C|9|‘, 


where C does not depend on b. For brevity we write 

II ■ li; = II ■ = (5) 

First, we state the result obtained in |Tren2] . He sets b = 

Theorem 2.1 { |Tre02j 1. Let conditions [H1-H3] hold. Then, there exist 
smooth functions 


and canonical coordinates {rj, h, r) such that the following conditions hold. 
• u = dg Adf + dh A dr. 
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• T] = ^+u{ri)T = cp+f, h = Hq- 

E{I)), where f denotes a function depending only on {I,p,q,e) satisfying 

/(j,o,o,o) = cg 


Define 

wl = h*- E{p*) - eU,{p*,^ + v{p*)T + 
For any a = {+, —} and {p* ,f, h* ,t) such that 

^-1^5/4| log^l ^ ^ ^^7/8 

|r| < c~^ c < \w^\ < c“\ 


( 6 ) 


(7) 


the separatrix map {p*,^*,h*,T*) = SA4{p,f,h,T) is defined implicitly as 
follows 


V* =V- ed^Q^{p*,f,T) -^^log 

r=e + /i"+ ed,.Q^ip*,f,T) + ^^\og 


A 


A 


h* =h- 


T* =T + t+ 


dh*Wo 

A 


log 

log 


A 


A 


+ O 2 
+ O 1 
+ O 2 
+ O 1 


( 8 ) 


a =a sgntfo, 

where \, and pr^ are functions of p*, i is an integer such that 

dh*Wo 


T + t + 


A 


log 


A 


< c 


-1 


(9) 


and Oi = log^e:, O 2 = log^e. 

Remark 2.2. Recall that the Hamiltonian has a normally hyperbolic invariant 
cylinder close to Aq = {p = q = 0}. Dynamics of restricted to A^ is 
often called inner dynamics. Dynamics of orbits belonging to the intersection of 
iy^(A£) and is called outer dynamics and thoroughly studied in jDdlLS0f)\. 

\DdlLS0Sl . 

It is helpful to have Figure]^ in mind. 

An heuristic explanation of the separatrix map is the following: every orbit 
starting in the fundamental region A and sufficiently close to the stable manifold 
fF^(A£) has three regimes: 

^One can show that / = ^ Hwq + e). 
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— approaching the cylinder A^; 

— evolution near 

— departing away from A^ and reaching the fundamental domain. 

In each of these three regimes we have: 

— Straightening invariant manifolds trivialized the first regime. 

— The Hamiltonians Wq approximate the evolution in the 2nd regime. 

— The splitting potentials describe the third regime with an error. 

It is reasonable to refer to the regime 2) as inner dynamics, since orbits near 
the cylinder A^ can be shadowed by orbits inside of the cylinder. 

It is reasonable to call outer dynamics to the regime l)+3). The outer dy¬ 
namics to the leading order is well described by the splitting potentials 0^^. 

Look now at the formula (7). The separatrix map has contribution from the 
Wq, which is the inner dynamics, and from 0^^, which are the splitting potentials. 


2.1 Formulas for auxuliary functions A, and 0^ 


The functions A > 0, > 0 and G M are dehned by the unperturbed 

Hamiltonian Hq as follows. Hypothesis H2 implies that both eigenvalues of the 
matrix 


A(J) 


-(9pgi7o (1,0,0) -(9qgi7o(l,0,0)\ 

<9ppi7o(l,0,0) <9pqi7o(l,0,0) ) 


( 10 ) 


are real and the trace of this matrix is equal to 0 for all I. We denote by \{I) 
the positive eigenvalue of this matrix. 

We denote by 7 ^(J, ■) : M —)■ the time parameterization of the upper 

and lower separatrices of Hq in the level of energy Ho{I,p,q) = We 

denote by a±{I) the left eigenvectors of the matrix A(J), that is, a+A = Aa+ and 
a_A = —Aa_, such that the 2x2 matrix with a± as columns has determinant 
equal to one. 

Then, we define 


/ -l-oo 

{-v{I) + diHo{I,-f^{I,t)))dt ( 11 ) 

■OO 

= \im [(a+(J),7^(/,-t))(a_(J),7^(/,t))e2^(^)*] . (12) 

t -,-- 1-00 

To dehne the functions 0”', we have to introduce some notation. We dehne 
the differential operator 

d = iy{I)d^ + dr. ( 13 ) 

Fix /d > 0. The inverse operator d~^ is defined on the “resonant” space 

Res = {/ : P X T”+^ : = 0 if \{k, u{I)) + ko\ < P/2} , 


13 



where are the Fourier coefficients of /. Then, for / G Res, we have that 
d~^f is well dehned and satishes 

As we have already mentioned, Treschev in |Tren2] chooses f3 = 

We dehne 


r (/, V?, r) = -d-^ 0, 0, r) - Hi(J, 99 , r)] . 


(14) 


where Hi is the Hamiltonian defined in (|^. 

Any solution of the unperturbed system lying on r^(/) can be written as 

= r^(/,(y9,r+ t) 


with 

r) = (/, v? + v{1)t + r), 7 ^( 1 , r)) 

where are solutions of 

lim x^(/,f)=0. 

t - > — 00 

By the definition of /i”" in ( [TT] ), 

= lim 

t - > + 00 

We define 

(p + ut + x‘^(/, ±cx)), 0 , 0 , f - r). 

Note that (p, r, t) tend to zero exponentially as t —> ± 00 . 

Then, 


‘d, -r) - ^ + /i'", -r) 


(p, r, t) dt 


-+00 ( 15 ) 

(p, r, t) dt. 


In |Tren2j these functions are called splitting potentials. The integral term is 
the classical Melnikov potential (also often called the Poincare function, see for 
instance IDGOOj h 

As we have explained the purpose of this paper is to refine the formulas given 
in Theorem 12.11 
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3 The finite harmonics setting 


To get more precise formulas for the separatrix map in the hnite harmonics set¬ 
ting, we consider different regions where it is dehned. These regions overlap each 
other. 

First £x some notation. Take a function / : T” x M"' x x T —> M with 
Fourier series 




Dehne Af as 

Af{f) = {ke : / V 0} 


and 

= {ke ■.k = ki + k2, ki, k2 e Afif)}. 

Consider the non-resonant region, which stays away from the resonances created 
by the harmonics in U 

Dehne 


NRg = {l:'ikeAf{H,)UAf^^'>{H,), \k-{u{I),l)\>(]}, (16) 

for a hxed parameter f3. The complement of the non-resonant zone is build up 
by the different resonant zones associated to the harmonics in J\f{Hi) 

Fix k G U then we dehne the resonant zone 

Re4 = {/:|4.M/),1)|<«. (17) 

The parameter (3 in both regions will be chosen diherently, so that the diherent 
zones overlap. 

We abuse notation and we redehne the norms in ([^ as 

II . Il* = II . ||(/3), c>F) = Ol = of. 

for hxed /3 > 0. Now we can give formulas for the separatrix map in both regions. 


3.1 The separatrix map in the non-resonant regime 

The main result of this section is Theorem |3.1| which gives rehned formulas for 
the separatrix map in the non-resonant zone (see (16)). To state it we need to 
dehne an auxiliary function w. This tc is a slight modihcation of the functions 
Wq given in (|^. 

Consider a function g{r],r). It is obtained in Section [tT by applying Moser’s 
normal form to Hq. This function g satishes g{r],r) = X{r])r where A is 
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the positive eigenvalue of the matrix (10) 
to the second variable for small r 


Therefore, g is invertible with respect 
Somewhat abusing notation, call g~^ the 


inverse of g with respect to the second variable Then, we dehne the function 
w by 

w{g*,h*)=g;\g*,h*-E{g*)). (18) 


Note that in the nonresonant regime, w does not depend on whether we are close 
to one of the unperturbed homoclinic loops or the other, as happened in Theorem 

EU 


Theorem 3.1. Fix j3 > 0 and 1 > a > 0. Let conditions hold 

for some s > 6, r > 8s + 2. Then for e sufficiently small there exist c > 0 
independent of e and a canonical coordinates {rj, h, r) such that in the 
non-resonant zone NR^ the following conditions hold: 

• the canonical form u = dg A df -\- dh A dr; 

• g = 1 + Ol{e) + f + v{g)T = ip + f^h = H^ + 0\{e) + 

O^^Hq — E{I)), where f denotes a function depending only on {I,p,q,e) 
and such that f = 0{w + e), f{I, 0,0, 0) = 0. 

• In these coordinates has the following form. For any a G , +} and 

{g*, h*) such that 


<\w{g\h*)\<ce, |r|<c-\ c < \w{g\h*)\e^^^*^^^ < c-\ 

the separatrix map {g*, , h*, r*) = SAis{g,^,h,T) is defined implicitly as 

follows 

g* = g — + C) 3 (e)|loge| 

r = ?+ di^'^{g,w{g* ,h*)) + dr,w{g* ,h*)[log\w{g* ,h*)\ -A d 2 ^''{g* ,w{g* ,h*))] 
+ Ol{e) I loge| 

h* = h- 0*^{e) 

t* = tA i A dhw {g*, h*) [log\w {g*, h*)\ A d 2 ^‘'{g*,w{g*,h*))] 

A Ol{s) I loge|, 

where w is the function defined in M* and are some func¬ 
tions and t is an integer satisfying ([^. The functions M* are evaluated at 
h*,T). 

This theorem is proven in Section 

^The subindex is to emphasize that the inverse is performed with respect to the variable r. 
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Remark 3.2. The change of coordinates in the above Theorem is e-close (in the 
C^-norm) to the system of coordinates obtained in Theorem 2.1 


The function M* and are defined in Lemmas 6.1 and 1.2 respectively, 
The functions are the generalizations of the functions /i”" and k®". Indeed, 
they satisfy 

+ 0;{r) and = k^{ 7]) + 0*{r). 

Moreover, the functions Mf’* satisfy 

= d^Q^ + 0*2 (w) , M['^ = dr0" + 0*2 (w), 


where 0'^ is the splitting potential given in (15) (compare with the formulas in 
Theorem 2.1). These last formulas are obtained in Lemma 1.3. 


Remark 3.3. We compare this Theorem with Theorem 2^ and the Remark af¬ 
terward. Notice that the inner dynamics in non-resonant zones can be made 
integrable and essentially decoupled from the outer dynamics. This is reflected 
in the fact that the rj and the h-components have no contribution from w as the 
inner dynamics is integrable. 

The splitting potential contributions and the integrable contributions to rota¬ 
tion of the angular components are essentially the same. 


Consider the following useful example: 

j2 p2 

He{I,T,P,Qfl) = y + y + + eHi{I,ip,t). 

Notice that the perturbation does not affect the pendulum. As the result the 
perturbed system is the product of the pendulum and a perturbed rotor. Notice 
that Hi(J,(p,f) = Hi{I,ip,t) and it is non-vanishing. Moreover, the splitting 
potential 0°^ vanishes. Since Wq, dehned in ([^, measures the trasition time i 
of the separatrix map SAie (see ([^), it has to be independent of angles of the 
rotor and time. This implies that in ([^ after substution of I,p,q,(p there is 
a cancellation of eHi with e-terms coming from (/ — rj) and {Hq — h). This 
cancellation is not so easy to see. 


3.2 The separatrix map in the resonant zones 


To give rehned formulas in the resonant zones Res^ (see (17)) we restrict to 
Hamiltonian systems of two and a half degrees of freedom. Namely, I and ip are 
one dimensional. 
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We consider the resonance (z/(/), l)-(/co, fci) = 0, (fco,/ci) E C 

where z/(/) is the freqnency map introdnced in ([^. We assnme that it is lo¬ 
cated at / = 0. Let A be the variable conjugate to time. Perform a change to 
slow-fast variables. Namely, consider the following change of coordinates 

(J, 9, D, t) = koip + kit, A - (19) 

which is symplectic. We obtain the following Hamiltonian 

H{J,9,p,q,t) = Hi [kol, ko^{9 - kit),p,q,t) . 

(kn All ^ — 

Note that in Res)^ ’ ^ and in slow-fast variables, the Hamiltonian Hi is time 
independent and thus Hi = Hi. It is defined as 


Hi(J,0) = 


2TTil9 




To state the next theorem, we recall the dehnitions of i/ in ([^ and of A, pA, tA 
and G®" in Section 2.1 We also consider a slight modification of the functions Wq 
in 

wo = h* -E{p*)-eHi{r]*,C)- ( 20 ) 


This definition is implicit since, as shown in Theorem |3.4| , depends itself on 
Wq. The function wq, as a function of p* and is independent of a. It certainly 
depends on a if considered as a function with respect to the original variables. 
As happened in Theorem |2.1 and unlike Theorem 3.1 the functions Wq depend 
on the homoclinic loop through the functions /i'^ (which appear in the definition 

ofD- 

We also define the functions 


R^’"(r/,e,<,r) = 


u{p) 


(Hi(r/,e + K^)r)-Hi(r7,0) 


Hv) 


Hiiv, ^ + h'"(^)) - Hi(r7, ^ + 


Kv) 


log 


K Wq 


A 


= Hiiv, ^ + h‘"(h)) - Hiiv, ^ + Hv)r)- 


( 21 ) 


Theorem 3.4. Fix /3 > 0 and 1 > a > 0. Let conditions [HI', H2,H3] hold for 
some s > 6, r > 8s-|-2. Then for e sufficiently small there exist c > 0 independent 
of e and a canonical coordinates {p, f, h, r) such that in the resonant zone 
Res^ the following conditions hold: 

• the canonical form u = dp A df + dh A dr; 
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• r] = I + Ol{e,Ho- E{I)), ^ + u{r]) = ^ + f, h = Ho + Ol{e, Ho - E{I)), 
where f denotes a C^~‘^ function depending only on (Lp,q,e) and such that 
f{I, 0, 0, 0) = 0 and f = O{wo + e). 

• In these coordinates SM^ has the following form. For any a G , +} and 
[r]*, h*) such that 

^ < ce, kl < c“\ c < |w;o(^*, Dl < c"\ 


the separatrix map {p*,^*,h*,T*) = SA4{p,^,h,T) is defined implicitly as 
follows 



h* 

=V 

+ 


) +eB^ 

yv,(, 

h,r) 



c 


V , 

+ ^log 

rfwQ 

A 

+ eB^ 


Wo,T,t) +0*{e) 


h* 

=h 

+ edrQ'’ 






T* 

=T -\-t 

+ ilog 

iHwq 

A 

+ eB^ 


U!o,T, 

t) +0*{e) 


a* 

=a sgnwQ 







where 








R« 

”(>) 

,^,h,T,t) = 


r)r + 


r) log 

kTWq 

A 

+ t)t 

R" 


, h, r, f) = 


t)t + 

92{v,^,h, 

r) log 

kTWq 

A 

+ 93{v,^,h,T, t)t 


for certain functions fi, Pi which satisfy fi,gi = 0*{1). Therefore, the 
functions satisfy 




_ C>*(log£). 


Recall that i/(0) = 0. Nevertheless, one can easily see that the function 
is well dehned even as z/ —)■ 0 since it has a well dehned limit. 

This theorem is proven in Section 


Remark 3.5. Vanishing the splitting potentials 0”' in the formula of the separa¬ 
trix map in the resonant zones one has the inner dynamics. Its non-integrability 
is given by the functions The splitting potentials G®" encode the outer dy¬ 

namics as in the non-resonant setting (see Remark 2.2). 
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We compare this Theorem with Theorem |2.1| and the remark afterward. No¬ 
tice that both the inner and the onter dynamics are nontrivial and £-conpled dne 
to resonant terms. Since transition time of the separatrix map is 0(\oge) this 
gives a term of order 0{e\oge) in the action components, which dominates the 
contribntion of the Melnikov fnnction. 

Note that in the aforementioned e-conpling between the inner and the onter 
dynamics vanishes in the generalized Arnold example and the Melnikov fnnction 
gives dominant contribntion (see Theorem A.2). 


3.3 Main steps of the proof and structure of the paper 

In this section we sketch the derivation of the separatrix map done by Treschev 
jTren2] . At each step we refer to the section where it is done in this paper. Recall 
that the separatrix map SM.^ is a retnrn map to a certain fnndamental region 
U“ U U’*' (see Fig. [^. Its compntation consists of six steps: 

• Moser’s normal form for the unperturbed pendnlnm near the separatrices 
(see the bine and the yellow part on Fignrej^ left and Section]^. 

• Moser’s normal form for the perturbed pendnlnm in the colored part of 
Fignre|^ left (see Section]^. 

• The transition map from one yellow region to another one in the variables 
given by Moser’s normal form of the nnpertnrbed system (see Section]^. 

• The regions U"*" (resp. U~) on Fignre|^ left are overlapping regions in the 
original coordinates (see Fignre|^ right). 

• Compnte the glning maps (Section]^. 

• Compntation of the composition of the transition map in Moser’s normal 
form variables with the glning map (Section]^. 


4 Normal forms 

4.1 Moser normal form close to the torus 

We start, as in |Tre02] . by performing the classical Moser normal form |Mos56] 
to the nnpertnrbed Hamiltonian Hq. To obtain a hnitely smooth version of this 
resnlt we apply a resnlt in |BdlLW96] . 

Lemma 4.1. Let Hq is satisfying [HI'] and mi. For I gV and {p, q) close 
to (0, 0) there exists a system of coordinates (/, s, x, y) such that 
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Figure 4: Gluing maps. 

• the change : {I,s,x,y) —> {I,(p,p,q) smooth and symplectic with 
i> (r-4)/5. 

• in the variables Hqo J^q = 'Ho{I,xy) is smooth. We write it as 

no{I, xy) = E{I) + g{I, xy) 
where g{I,xy) = X{I)xy + 02{xy). 

Proof. In Lemma l. |Tre02] the assumption is that Ho is analytic in p, q. To relax 
this assumption we use Theorem 1.1 jBdlLW96j. Recall that bv Remark 0 we 
can assume that the saddle is at (p, q) = (0,0) for all / at the expense of loosing 
two derivatives. 

Let / be a diffeomorphism with the hxed point at the origin /(O) = 0. 
Let iV be a symplectic polynomial map such that A^(0) = 0 and the of / 
and N coincide at 0, i.e. f{0) = D^N{0), j = 0,..., fc. Let 1 < i < kA — B 
and r — 2 > + 4 for some integer £. Then N and / are conjugate, i.e. there 

is a diffeomorphism h such that h~^ o f o h = N near the origin, h(0) = 0. 

In our case for each hxed I we have a 2-dimensional symplectic map / with a 
saddle hxed point. By the remark after Theorem 1.1 of |BdlLW9^ . R = 1/4, B = 
1-2R = 1/2. Thus, 1 <£< ^ < ^. □ 

We consider the expansion in e of the perturbation of the Hamiltonian ([^, 
namely. 

Hi = Hii eHi2 0(8“^) 

and dehne 

n, = no + eHi + e^n2 + O(e^) (22) 

with 

'Ri(/, s,x,p,f) = Hiio{Ho{I,s,x,y)A), %(/, s, x, p, f) = Hi 2 o{HoiI,s,x,y)A)- 
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Since we assume that r/8 — 5/4>swe have that 1-L^ is Hi is and H 2 is 
Moreover O(e^) = Ocs- 2 {e^). 


4.2 Normal forms near the separatrices 

We extend the Moser normal form to the region where \xy\ is small (see the 
shaded/colored region Fig. |^on the right/left respectively). This normal form 
applies to both the non-resonant and resonant regimes. Since we want to make a 
second order analysis of the separatrix map we need a more precise normal form 
than in |Tren2] . In this normal form there are two sources of error. 

• Expansion in the small parameter e: we need to perform two steps of normal 
form instead of one to reduce the size of the remainders. 

• Powers of xy. Treschev only performs normal form to remove the terms in 
the perturbation which are independent of the product xy. He takes 

loge| < \xy\ < 

We want to remove terms up to the first order in xy. Assume that 

^ \xy\ < £ with 1 > a > 0. (23) 


We perform the change of coordinates by the Lie Method. First, we proceed 
formally and then we compute the estimates. We consider the expansion of the 
Hamiltonian PL given in (22) and of a Hamiltonian of the form eW = eWq + e^Wi. 
Call $ the time-one map associated to the flow of eW. This change is symplectic. 
Moreover, 


H, o $ =H, + e{H, W} + £^{{n, W}, W} + 0 (e^) 

=Ho + £(Hi + {Ho,IFo}) 

+ ({Ho, hhi} + (Hi, Wo} + {{Ho, Wo}, Wo} + H2) + O {e^) . 


Now we look for suitable Wo and Wi. 

First, compute Wq and then Wi. To compute Wo, split the Hamiltonian Hi, 
defined in ( [22| ), in the following way 

Pii{I,s,x,y,t) =Hi(/, s,f) +PL^^\l,s,t) +PL^‘^\l,s,y,t) + PL^^\l, s,x,t) 

+xy (H 2 (/, s, t) + PL^^\l, s, t) + PL^^\l, s, y, t) + PL^^\l, s, y, t)) 
+o;{xy) 


where 


ii,{i,s,t)= ^ 




a,X(^.0,0)e 


27 ri/c-(s,t) 


(24) 
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Ip is the bump function introduced before Q and 

s, t) = 5, 0 , 0 , t) - Hi(/, s, t) 

s, X, 0, f) — s, 0, 0,f) 
Uf\l,s,t) = 7 ^ 1 ( 1 , s, 0 , 1 /, f) - 'Hi(/, s,0,0,f) 
s, t) = d^yH,{I, s, 0,0, t) - H2(/, s, t) 
Uf\l,s,t) = d^yHi{I,s,x,Q,t) - 5a;j/'Hi(/,s,0,0,f) 
'Hf\l,s,t) = da:y'Hi{I,s,0,y,t) - d:^y'Hi{I,s,0,0,t), 


where Hi is dehned in (4). The functions Hi, 'H^i \ j = 1, 2, 3 are whereas the 
functions 'Hi \ j = 4, 5, 6 are 

The next lemma contains the hrst step of the normal form. In the next two 
lemmas we denote by {•, ■}(x,y), the Poisson bracket with respect to the conjugate 
variables {x,y). 


Lemma 4.2. There exists a smooth solution Wo{I, s,x,y,t) of the equation 

3 6 

(z/(/) + dig{I, xy)) dsWo + d^Wo + {//(/, xy), ITo}(.,,) + U? +xyY, = 0- 

j=i i=4 

The functions Wn satisfies 

Wo = o*{r^), 

where O* is defined in &■ 

Proof. We take Wq = and solve the equations 

+ dig{I, xy)) dsWl + dtWl + {//(/, xy), = 0, j = 1, 2, 3 

(z/(/) + dig{I, xy)) + dtWi + {g{I, xy), + xyUf = 0, j = 4, 5, 6. 

Each equation is solved as follows. For the hrst and fourth ones, we just have 

(z/(J) + dig{I, xy)) + dtW] + -hS') = 0 
{v{I) + dig{I, xy)) + dtW^ + xyH^^ = 0 . 

Thus, we invert the operator 

d := {iy{I) + dig{I, xy))ds + dt (26) 
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by using the Fourier expansion and inverting for each Fourier coefficient. Note 
that this operator and the operator 9 in (13) satisfy d — d = 0{xy). Moreover, 
recall that 'Hq is and therefore so is g. Then, Wq'"^ are and satisfy 

and = 0*(xy/3-^). 

For the others, we use the characteristics method to obtain 

rO 


^0 = 




’ —OO 

/*+oo 


s + (z/(J) + dig{I, xy))t', ^ 

nf\y, s + (i/(J) + dig{I,xy))t\xe-^^3{i,xy)t' 

Jo 

pO 

W! = -xy / V!i\l,s^{u{l) + dig{l,xy))t\ye^^3{i,xy)t' 

J —OO 
r+oo 

Wq =xy s + {v{I) + dig{I, xy))t\ , t + t') dt'. 


Thus, they are all 

Now we solve the second order equation. Dehne 

^2 = {Hi, IFo} + {{Ho, IFo}, PFo} + H2, 


□ 


where 7 ^2 is Hamiltonian dehned in (22). Using the equation for Hq, given in 
one has that H 2 is C^~^. 


4.2 


Lemma 

Let H 2 (/, 0, 0), A; G denote the Fourier coefficients of H 2 in s and t. We 
split H 2 in several terms, as done for 'H}, in the following way 

n 2 {I,s,x,y,t) =H 3 (J,s,t) +'Hf\l,s,y,t) + "Hf ^(J, s, x, f) 

+ 0*{xy) 


with 


and 


H3(/,S,t)= V' 


\ ^k/T n n\„ 27rik-{s,t) 


f3 


H^(/,0,0)e^ 


(27) 


H^^^(/, s, t) = 772 ( 7 , s, 0, 0, t) - H 3 (/, s, t) 

H2^^(/, S, t) = H2(/, S, X, 0, t) - H2(/, s, 0,0, t) 

s, t) = '^ 2 ( 7 , S, 0, y, t) - H2(7, s, 0, 0, t). 

All these terms are C^~^. 
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Lemma 4.3. There exists a ^ smooth solution Wi{I, s,x,y,t) of the equation 


(z/(/) + dig{I, xy)) dsW^ + dm + {9{I: xy), = 0. 


i=i 


Moreover, Wi = 0*{(5 ^), where O* is defined in 0. 

we take Wi = solve the 


Proof. As in the proof of Lemma 
equations 


4.2 


{u{I) + djg{I,xy))dmi + dtWl + {g{I ,xy),Wl}^,,,ym = 0, J = 1,2,3. 

Each equation is solved in the same way as the hrst order in Lemma |4.2 One 
can see that 'H 2 satishes 7^2 = 0*{/3~‘^). This implies that Wi = 0*{l3~^) (one 
can have more precise bounds for each Wl). □ 


Denote by 

$ : (/,s,x,^,f) I—^ {I,s,x,y,t) 

the time-one map associated to the flow of the Hamiltonian eW = sWq + e^Wi. 
This change is and symplectic. In the next two lemmas we analyze the 
change of coordinates and the transformed Hamiltonian. 


Lemma 4.4. The change <h zs C® and satisfies the equations 

1=1 + eMl + e‘^M^ + O* (e^/3-^) 

= 1 + sdmo + {dsWi + {dsWo, hho}) + O* 

s=s + eMm + O* 

= s- ediWo - {diWi + {djWo, Wq}) + O* 
x=x + eMf + e^Mf + O* 

= x + edyWo + {dyWi + {dyWo, ITo}) + O* 
y =y + eMf + e^Mf + O* 

= y- edmo - {d^i + {d^o, hho}) - O* 


Moreover, 


= 0*(/3-‘) 

, z = I, x,y 

and Ml 

and 

Mf 


z = I, x,y 

and Ml 

We also have 






xy = xy + eM[ -f + O* , 
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with 


Ml = and Ml = 0*{(3-^). 

The inverse change is of the same form, that is 

T=I- eMi + + O* (e^/3-^) 

s=s- eMl - e^Ml + O* 
x=x- eMf + e^Mf + O* 
y=y- eMf + e^Mf + O* 

and ^ 

xy = xy- eMl + e^M^ + 0*{r^e^). 

The terms M| satisfy the same estimates as M|. 

Proof. It is enough to recall that 

1 = 1 + e{I, W} + £'{{/, W},W} + .... 

To compute the remainder, one has to estimate 

{{{z, ITo}, hho}, ITo} + hho}, ITi} + {{^, ITi}, ITo}, z = I,s, x, y. 

Using the estimates for hUo and lUi, one obtains the bounds for the remainder. □ 

Now we can apply this symplectic change of coordinates to the Hamiltonian 
PLe given by Lemma |4.1 

Lemma 4.5. The Hamiltonian Pie o ^ is and is of the following form 

Pie o ^{T,s,x,y,t) =E{T) +g{T,x]f) + eHi{T,s,t) + exyH 2 {T,s,t) 

+ H 3 (J, 's,t) + 0* + e‘^l3~‘^xy + e(xy)‘^) . 


5 Transition near the singularity dynamics in 
the normal form 


We compute the equation associated to the Hamiltonian given in Lemma |4i5| We 
drop the hats to simplify notation. Following |Tre02] . consider a region for the 
initial conditions of the form 

17* = {{I,s,x,y) : c* < |x| < cL\ Cq + \xy\f log^{xy) < \xy\ < At*} (28) 

and a final time t with 
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Take, as in |Tren2j . c* and cq independent of £ and k* ~ 

The Hamiltonian obtained in Lemma l4^ has different formulas in non-resonant 
and resonant zones. We hrst analyze it in the non-resonant zone and later in the 
resonant one, which are defined in (16) and respectively. 


5.1 The non-resonant regime 


Recall that we study trigonometric perturbations. In the non-resonant zone (16), 
analyzing (j^, (24), and (27) we have that Hj = 0 for j = 1,2,3. Therefore, we 
have the Hamiltonian 

Tie o ^(T,s,x,y,t) = E{T) +g{T,^) + O* , (29) 

which is and is integrable up to order 3. 


Lemma 5.1. Suppose that for some (R, s*, a:*, ?/*) G (see (28 )j and t G M, 


where p = \x^y^\. 

Then, 

s{t) = s* + + djg{r,p)) t + O* -1- s^j3~^p + ep^) log^ p 

lit) = r + O* (£3/3-4 + I logpi 

x{t) = X,e-^3{I*,p)t + Q* (^3^-4 ^ ^ ^^ 2 ^ 

yit) = (1 + O* (£3/3-4 + £2^-2p + ep^)) 

x{i)y{t) = x*y* + O* + s^l3~‘^p + ep^') \ log p\. 


Proof. The proof of this lemma is a direct consequence of the particular form of 
the equations associated to Hamiltonian (29). Indeed, one can easily see that 


d 

dt 


(xy) = O* (£3/3 ^ -1- £^/3 ‘^xy + eixyY) . 


Therefore, one can easily see that 


xit)yit) — x*y* I = O* (£3/3 ^ -\- £^/3 ep^) \ log p\. 


Taking this into account, we have 

i = 0* (£3/3-^ + e^l3-^p + £p2) 

which leads to the formula for I(t). Using that I is almost constant, one can 
easily deduce the formulas for the other variables. □ 
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5.2 The resonant regime 

To analyze the resonant regime recall that we focus on the case of two and a half 
degrees of freedom. Namely, s G T and / G M. We perform a change to slow-fast 
variables. This leads to a Hamiltonian which is almost a first integral, namely, 
its time dependent terms are small. 

Fix {ko, ki) G I?. Assume that the resonance (zz(/), 1) ■ (/cq, ki) = 0, (/cq, ki) G 
is located at / = 0. Call A the variable conjugate to time. Then, 

the change 

(J, 0, da) = ( 7^, kos + kit, A - p-I, t 
\ko ko 

is symplectic. Applying this change, one obtains the following Hamiltonian. We 

drop the hats to simplify notations. 

V-iJ, 0, X, y, t) = E{J)+7j{J, xy)+el-Li{J, 6, xy) + 0* -I- e^l3~‘^xy + e^xy)"^) . 

where 

E{.J) = E{koJ) + kiJ, 

satishes ^ 

djE{0) = 0, g{J, xy) = g{koJ, xy) 

and 

UiiJA^xy) = 

N 

E fTliAoJki) , , YliAoJki) , , ^(jkojki) , 2-Kiie 

{koJ)+xyU2 {koJ)+eU^ {koJ)je W 

j=-N 

We use this system of coordinates to analyze the flow in the resonant zones. 
Recall that by construction z/(0) = 0. Since we consider /3 > 0 hxed and in order 
to avoid cluttering the notation, from now on, we do not keep track of the fd 
dependence of each estimate. We also assume p < 

Lemma 5.2. Suppose that for some {E,9^,x^,yA G and t G M, 

where p = \x^y^\. 

Then, 

e{i) =e* + (HJ*) + dig{J\p)) t + eEi{.J\ e\t) + O* 

J(t) =.r + eGi{J\ e\t) + O* 

x(i) + O* 

y(i) =y^e^ 2 ff(J*,p)t+e'^(J*r,t) ('i + Q* 

=x*y* + O* , 
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where 


N 


' k=-N,k^O 
N 


21719 


djH,{J)t+ Y. 


2'Kikh'{J) 


^2niku{J)t _ 

2'Kikv{J) 




(30) 


k=-N,ki^0 

and 

$(J,0,t)= [ {U 2 {J,e + iy{J)t) + drjg{p,J)Gi{J,e,t))dt. 
Jo 

Thus, G has zero average with respect to 9 and F and G satisfy 
F = 0*{\oge) and G = 0*{1). 


(31) 


Proof. From Lemma 4^, we have the equations 

9 = u{J) + djg{J,xy)+edjH,{J,9)+exydjH2{J,9)+e^djH3{J,9) 

+ O* + e‘^/3~‘^xy + e{xy)'^) . 

j = - edeHi{J, 9) - exyden2{J, 9) - e'^deHsiJ, 9) 

+ O* + e‘^f3~‘^xy + e{xy)‘^) 

y = drg{J,xy)y + eyTl2{J,9) + Ol{el3~^ ,xy) 

X = — drg{J, xy)x + exli 2 {J, 9) + 02{e/3~^, xy). 

We also have 

j^{xy) = O* + e^/3~^xy + e{xyf) . 

Now we compute estimates for this flow. Call {I*, s*,q*,p*) to the initial point. 
Recall that p = x*y*, that we omit the dependence on fj and that we assume 
p<e. Then, 

\xy — x*y* I = O* log e) . (32) 

This implies that the hrst orders of 9 and J are independent of x and y, since 
only depend on p. Indeed, we have 

9 =i^{J) + djg{p, J) + edjlii{J, 9) + 0*{e'^) 

j = -edeHi{J,9) + 0*{e'^). 

These equations can be solved perturbatively in powers of £. We look formally 
for solutions of the following form and we take t* = 0 

9{t) = 9* + iy{r)t + djg{r, p)t + eFi{r, 9*, t) +h.o.t 
J{t) = J*+ eGi{J% 9*, t) + h.o.t 
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Plugging these expressions into the equation and recalling that g{J*, p) = djg{J*, p) 
we obtain the following equations 

Fi = iy'{r)Gi + djU^r, e* + p{r)t) 

Gi = -deRi{J\e*+ p{.r)t) 

We can hrst solve the second equation, taking 

Gi(J*,r,t) = - [ deUiiJ^e* + iy{J*)t)dt, 

Jo 

which leads to the formula of Gi stated in the lemma. 

Plugging the expression of Gi into the Fi equation , we obtain 

Fi{j*,e*,i) = -F{j*) [ Gi{r,e*,t)dt 

Jo 

+ [ djHi{r,e* + u{r)t)dt. 

Jo 

Now, expanding every term in Fourier series and integrating we get the formula 
for Fi. From their dehnition we can deduce that the functions Fi and Gi satisfy 


F, = 0*{\oge), Gi = 0*{l). 


Analyzing the remainders, one can easily compute the size of the higher order 


terms in the evolution of J and 9 given in Lemma 5.2 


Now we analyze the flow for the x and y variables. Recall that the conditions 
on the initial y^, and on the time t imply 


G ^ ^d2g{J*,p)t ^ 

P ~ ~ pet 


(33) 


Using the almost conservation of xy and the formulas for 9{t) and J(t), we have 


y = {drg{p, J{t)) + eH 2 (J(t), 9{t))) y + 0i(f) 

X = - {drg{p,J{t)) + etl2{J,9{t))) x + 02(t). 

where 0* = 0*{e^\oge). 

Thus, 

J a 

"‘The terms Fi does not appear in |Tre02) 
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Using (33), 


x{i) = x^e-^o{drg{p,J{t))+eU 2 {J{t),e{t),t))dt ^ C>*(^ 2 log 2 ^)_ 

Now we expand J{t) and 6{t) with the formulas already obtained. Then, 
x{i) = 

with 

= r (H2{J%e* + iy{J*)t) + drjg{p,J*)Gi{J%e%t)) dt. 


Using this estimate and (32), we can deduce the following formulas for j/(f), 
y(i) = y^e^r9{p,r)i+eHJ*,e*,i) ^ c>*(eiV6)) _ 


□ 


6 The gluing maps 


We compute the gluing maps. First in Section 6T, we consider the unperturbed 
Hamiltonian "Ho given in Lemma [4.1[ Then, we consider the full Hamiltonian in 
the non-resonant setting (Section |6.2[) and resonance setting (Section |6.3D. 


6.1 Unperturbed gluing maps 

The unperturbed gluing maps are computed in Section 5 of jTren2] . Consider the 
gluing maps 


5'*' is small,a; > 0} —> {y > 0, |x| is small} 

S~ is small,a; > 0} —> {y < 0, |x| is small}. 


(34) 


The gluing maps of the unperturbed system must satisfy the following properties 

(i) are symplectic 

(ii) preserve I and xy. 


Lemma 6.1 f |Tre02j L Properties (i), (ii) and (34) imply that the gluing maps 

(U, s^, x^, y^) = S^{I, s, X, y) 
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are of the form 


1^ = 1 

= s + xy) 

= x-l^-dr^^(hxy) 

for some functions . 

In |Tre02j it is also shown that di^^{I,r) = jx^{I) + 0{r) and = 

K^{I) + 0{r). Notice that the maps are C*. 

Moreover, since p = xy^ x^ = 0{p) are always small. One can consider also 
the inverse change for {x,y). Since xy = it is given by 




So, y also is small since y = 0{p). 


6.2 Gluing maps in non-resonant zones 

Now we express these gluing maps in the normal form coordinates (J, s', x, y). We 
use the formulas for the normal form coordinates obtained in Section 14.21 

The gluing maps have the following form. Each term can be expressed in terms 
of the Hamiltonian W associated to the normal form variables from Section 14.21 

=I + e{Mi oS^ - Ml) + e^{Ml o - Ml) + 0*{e^(5-^) 

^ =s + + e{Ml oS^- M/), x^+ £(M[ o 5^ - M[)) 

+ e{Ml oS^ - Ml) + 0*{e^/3-^) 

y^ =x-^e-^rMiT,xy) ^ 

j (x^+ £(M[ o 5^ - M[)) + 0*{p^ + pe(5-^ + e^(5-^) 

Moreover, 

y^x^ = yx + e{Ml o - M{) + e^{Ml o - M^) + 0*{e^(5-^). 


6.3 Gluing maps in resonances 

We consider the slow-fast variables {9,J,x,y). In these variables, we have the 
gluing map 

= S^{J,6,x,y) 
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with 


J^ = J 

6*^ = 6* + X?/) 

a;± = 


where 

$^(/,r) = ^^{koI,r). 

From now on, we drop the tilde in to simplify the notation. We also abuse 
notation and we consider the normal form change of coordinates given by the gen¬ 
erating function W expressed in slow-fast variables (recall that all these changes 
of coordinates are symplectic). 

Now we express the gluing map in the normal form (slow-fast) coordinates. 
As before, for the hrst four coordinates we have 

J± =J + £(M/ o - M/) + 0*i€^) 

=e + dj^^{J + e{M( oS^ - M(),xy + e{Ml o - M{)) 

+ e{Ml oS^ - Mf) + 0*{e^) 

+ C>*(e) 

-± £( m [ o - Ml)) + 0*{e^) 

and ^ 

y^x^ = yx + e{Ml oS^- Ml) + e\Ml o - M^) + 0;{s). 




7 The separatrix map 

We use the results in Sections |5.1| an d 6.2 to look for formulas of the separatrix 
map in the non-resonant regimefl6| . First we compose the flow in normal form 
coordinates and the gluing map. We obtain the separatrix map in the (/, s, y, x) 
coordinates. Later we look for a good system of coordinates which will transform 
{y, x) to a certain symplectic flow-box coord inat es around the former separatrix. 


Consider the flow gl analyzed in Lemma 5.1 and the gluing map analyzed in 


Section 6.2 We have the following 

Lemma 7.1. The composition of the two maps = g^o Gf is given by 


^{I,s,y,x) = (J'/(/,s,?/,x), J's(/,s,|/,x), J'j/(/,s, 2 /,x), j;(/,s, 2 /,x)) 
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with 


7I =I + eV( + + e^/3-^p + ep‘^)\\ogp\ 

7s =s + di^^{7xy) + (z/(/) + dig{I,xy))t + eVl + O* {e^ l3~^)\\og p\ 
7y =x~^ eyip{—drg{I + eV[^xy + eV{)i — dr^^{I^xy))x 
{l + Ol{e,p)\\ogp\) + Ol{e,p) 

7x =exp{drg{I + eVl,xy + sVl)i + dr^'^{I,xy)))x 

x{xy + eVl) (l + C> 2 (£,p)|logp|), 


where 

V[ =Mi oS^ - Mi 
Vi =Mi oS^ -Mi 

and 

VI = Ml oS^- Ml + xy) {Mi o - Mi) + 

+ xy) {Ml o - Ml) + tdi{u{I) + dig{I, xy)) {Mi o - M/) 

+ tdlig{I,xy) {MloS^-Ml) 

Moreover, 

7,s7y = xy + eVl + e^Vl + O* {el3-\ p) , 

where 


VI =Ml oS-Ml 
V; =Ml oS-Ml 


Then, 


Vi,Vl = 0*{r^), Vi,Vl = 0*{r^), and VI = 0*{l)\\ogp\. 


The proof of this lemma is a direct consequence of the results in Lemma 5.1 
and Section W72\ 

We look for a coordinate change near the former separatrix such that formulas 
for the separatrix map are as simple as possible. In comparison with |Tre02j . we 
want to point out two main differences. First, since we are away from resonances, 
we have Hj = 0, j = 1,2, 3. This simplihes the formulas. On the other hand, we 
want to have a more precise dependence on xy since we are doing a higher order 
analysis. This second fact implies that we have to slightly modify the change of 
coordinates. 
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We look for a symplectic change 


(s\ 


n\ 

I 

= T 

g 

y 


T 

\x) 


\h) 


The function I i—)■ g{I,r) introduced in Lemma 4.1 is invertible with respect to 
r in a neighborhood of r = 0 (recall that drg{I,0) = A(/) > 0, see Lemma 4.1). 
We have denoted by g~^{I,r) the inverse function with respect to the second 
coordinate r. We consider the following generating function. 


S{v,s,h,y) =ris + g/{ri,h- E{ri))\og\y\. (35) 


This generating function induces the following change of coordinates 

/e\ 

V 

T 

\hj 


(s- drg^ g{I, xy)){u{I) + dig{I, xy)) log \y\\ 

I 

log hi 


V 


drg{I,xy) 

E{I)+g{I,xy) 


and 


(s\ 


(i\ 


I 

= T 

g 

_ 

y 


T 


\x) 


\h) 



a exp( 


g 

_T_' 

drgr^(r],h-E{r ]))' 


(36) 


drgr^{rj,h-E{r])) ^/ 

Note that the y component does not get modified by this change of coordinates 
and that this change of coordinates does not depend on £. 

We express the separatrix map in these coordinates. We use the change T 
to write down the formulas. From now on we omit the dependence on (3. Recall 
that /3 > 0 is a hxed parameter independent of £. Due to the previous analysis, 


smoothness of the separatrix map obeys the estimate in Theorem 3.1 
Lemma 7.2. The separatrix map, has the following form 
g* =g + eV{ o T + e^Vi o T + 0;(e, p) \ log p\ 


C + di^^ig, g^ \g,h- E{g))) - 


iy{g) + djg{g, g^ ^{g,h- E{g))) 


drg{v,9T^{v,h- E{g)) 

+ C>i(^ + p)(|log^l + |logp|) 

h* =h + e{[E{g) + dig{g,g;\g,h-E{g)))] Vi oT 

+drg{g. g;\g, h - E{g)))V[ O T) + e^v!f + 0*(£, p) 

+ * + dMri, 9 Al,h-E{ri)) « + O;). + p)(| log.| + | logp|). 


logi? 
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where 


^ _ ^dr<s>^{wgr \v,h-E{v)) h* - E{7]*)) + O;) 


(37) 


= i^iv) + digih^gr \g.h- e{ 7])))) vi o t 

+ drgiri, g-^{ri, h - E{vi)))Vl o T 
+ ^ [dii^{g) + djg{g.gr^{g.h - E{ri)))) {vi oT)^ 
+ \dlg{7], g;^{g, h - E{g))) (P[ o T)^ 

Proof. We start with the g component. We have that 


V* = y* = y + £Pi + + Osie, p) \ log p\ 


Then, it is enough to apply the change T defined in (36) to obtain the formula 
for g*. 

For the h component we have, 


h* = E{r)+g{I\x*y*) 

= E{I)+g{I,xy) 

+ (i/(/) + dig{I, xy)) {eVi + e^Vi) + drg{I, xy) {eVl + e'^Vl) 

+ y + d]g{I,xy)) {ViY + ^—dlg{I,xy) {Vlf + 0^{e,p)\ logp| 

Apply the change T, one obtains the formula for h*. 

To compute the r component we use the following identity, 

g;\g% h* - E{g*)) = g;\g, h - E{g)) + eVl + 


We also have 
* log \y* 


drg{I*, x*y*) 

log \y\ + drg{I + eV(,xy + eVl)i + log (xy + eVl' 


drg{I + eVi,xy + Vl) 


+ C>2| logp|. 


Then, we obtain 
T* = T + i + 


drg{v,gr^iv,h- E{g)) 


logi? + C>;;(£ + p)(|loge| + |logp|). 


where R is the function introduced in (37). 
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Proceeding analogously, one can compute the ^ component, 


C =^ + di^^iv,9r\v,h- E{t]))) - 
+ Ol{e + p){\\oge\ + |logp|). 


iy{r]) + digjr], \v,h- E{r]))) 
drg{v,gr^{v,h- E{r])) 


\ogR 


Note that dig{ri,g — E{g))) satisfes 


dig{v,gr\v,h- E{g))) = Ol{e,p). 


Therefore, 


r g^ \v,h- E{p))) - 

+ C>t(£ + p)(|log£| + |logp|). 


_Kh)_ 

drg{v,gr^{v,h- E{p)) 


\ogR 


This completes the derivation of the separatrix map in the non-resonant case. □ 


Theorem |3.1 is a direct consequence of this lemma. We dehne the hrst orders 

as 


in the action components in Theorem 3.1 


= VloT 

Mf’" = [z/(p) + djg{p, g;\g, h - E{p)))] VioT + d^v, 0)iP[ o T. 


The second orders can be also easily dehned. 

We hnish this section by identifying the order £ of the separatrix map in terms 
of the Melnikov potential. This allows to compare our results with the results in 

|Tren2] . 

Lemma 7.3. The functions and satisfy 

+ 0(h - E(p)), + 0(h - E(p)) 


where 0°^ are the splitting potentials defined in (15). 


Proof. We explain how to prove the statement for the p component. The other 
one can be obtained analogously. 

In Lemma 15.11 we have seen that 


= V(oT = M(oSoT-MloT. 

We start by analyzing the function M( = dgW^ given in Lemma 4.4 
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By the definition of the fnnction Wq in Lemma |4.2[ we have that Wq = 
Woo + 0{.xy) where Woo is defined as follows. Recall that dig {I, r) = 0{r) and 
drg{I,r) = A(/) + 0{r) with A(/) > 0. The fnnction ILoo can be split as 


with 


IToo = f^oo + W^oo + ^oo> 


IToo = do^-H? 


= - 


' —oo 
/‘ + 00 


W^oo = - 


s + + t') dt' 

s + v{I)t\ + t') dt' 


where "Hj are the Hamiltonians introdnced in (25) and 


d-\f) = d 


-1 


E 


yfcAo(j)g: 


27Ti{ks-\-kot) I _ 


(k,ko)£Z' 


n + l 


E 




2711 {kiy{I) + ko) 


^27ri{ks-\-kot) 


(fc,fco)eZ^+i 

which is well defined in the non-resonant zone. These functions are the first order 
in xy of the funtions Wq, j = 1, 2,3, introduced in the proof of Lemma 


4.2 


we 


Therefore, M( = dsWQo + dsWQo + dsWQo -|- 0{xy). Now, in Section 6.1 
have seen that y^,x = 0{xy). Moreover, from the definition of the Hamiltonians 


in (25) one has that Til = 0{x) and = 0{y). Thus, 

M(oS- Mi = dMlo ° - dsWl^ - dM^o + dsW^^ o S + 0{xy). 


Moreover, as shown in |Tre02j and recalled in Section 6.1, the functions satisfy 
(?/$=*"(/, r) = + 0{r) and = k^(/) + 0{r). Then, we obtain 

MioS- Mi = dsWioil, s + /i"(/), t) - dsW^oik t) 

+ dMiioik s + /i"(/), y, t) - dsW^oil. s, 2 /, t) + 0{xy) 

= dsWUk s + /i"(/), t) - dMUk s, t) 

+ dM!o{k s + (/^"(/)2/)-\ t) - dsWl{I, s, 2 /, t) + 0{xy). 

Now it only remains to apply the change of coordinates T. The change T satisfies 

h, K r) = (^ + iy{g)T, g, 0) + 0{h - E{g)). 
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We apply this change of coordinates to each term. We obtain hrst 


s + 0) - d.WUh *>.•))I= 

(, + v(ri)T, 0) - { + v(:ri)T + 0), 


where d is the fnnction introdnced in (14). 

For the two other terms, it is enongh to use to see 

0) 1 {I,s,y, 0 )=r{ri,i,T) ~ 

dsWi^oiv, ^+i^iv)r + 0 ). 


Lemma 7.4. The following identity is satisfied. 


^{V, ^ + ^{v)^ 0) - d{r], ^ + u{r])T + p'"(77)r, 0) 

+9sW^Q{y, ^ + iy{T])T + ia^{T]), ^(r;)-0) = 

^{v, C, -r) - ^{v, ^ 

+dsW^Q{7],^,a, -r) + dsW^Q{y,^ +fi^{7]),K{7])-^a,-T). 


Proof. This lemma follows from the definition (14). One can expand both sides 
into Fourier series and match them. □ 

□ 


From this lemma, one can easily deduce the statement of Lemma 7.3 


8 The separatrix map in the neighborhood of 
resonances 

In the resonant regime we only compute the system up to first order. Thus, we 
follow closely |Tre02j . The main difference is that our resonant region is much 
larger than in |Tren2j . Indeed, our fi is hxed indepedent of e, whereas in |Tren2] 
he considers fi = Recall that we are not keeping track of the dependence on 
fi and that we have assumed p ^ e. 

First, we compose the flow in the normal form coordinates and the gluing 
map. We obtain the separatrix map in the {J,9,x,y,t) coordinates. 

We denote the composition of the two maps o by 

J^{J,e,x,y) = {J^j{J,e,x,y),J^e{,J,9,x,y),J^:^{J,e,x,y),J^y{J,e,x,y)). 
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Recall that the map is independent of t because we choose the initial time as 
t = 0. 

For the J component one can easily see that 

J^j = j + eV( + 


where 

V^{J,9,x,y) = Gi{J,9 + + M/ oS{J,9,x,y) - M^{J,9,x,y) 

where the function Gi has been introduced in Section [RS and the function has 
been used to define the gluing maps in Section 6.3 From Lemma 4.4, Lemma 


5.2 and the dehnition of the gluing map in Section 6.3, we can deduce that 


V( = 0*{l). 

For the 9 variable we have 


70 =9 + dj^^{J,xy) + (z/(J) + djg{J,xy))t 

+ eVl + 


where 

Vl =Fi + Mf o 5 - Ml + (J, xy) {M( oS- M() 

+i {dju{J) + djgiJ, xy))) {M( oS - M() + td%g{J, xy)){Ml oS - M[). 


Using Lemmas 4.4 and 5.2, one can check that Vl 
For the x and y components we have, 


0*(l0g£). 


7x = y ^ exp(— {drg{J + sV(, xy + eVl)t + <F(J, 9, t) + xy)))x 

{i + o;is)\\ogs\) + o;{e) 

7y = exp((9r5'(</ + xy + eV{)t + <F(/, 9, i) + xy)) x 

y{xy + eVl) (1 + C>2(£)| loge|) + C>3(£). 


Moreover, 


where 


7x7y = xy + eVl + 0{e^ 


V[ = Ml oS 

by Lemma 4.4, satishes VI = C>(/3“^). 


M[, 


To have simpler formulas for the separatrix map we consider a change of 
coordinates slightly different from the ones used in |Tre02j and in Section The 
reasons are the following. On the one hand, we have to deal with a rather large 
resonant region since it has width of order one with respect to e instead of order 
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as in |Tren2] . On the other hand, we do not need to consider higher orders 
as in Section since in the resonant regime we just analyze the hrst order of the 
separatrix map. As in the single resonant regime, we assume that \xy\ < e. We 
consider the change defined by the generating function 

S{J,i,x,T) = Ji +E{J)T + axe^^^'^^+ e ! Hi(J, ^ J)s)ds. (38) 

Jo 

Note that in the core of the resonance i^{J) ~ 0 and, therefore. 


Hi(J,e + z/(J)s)ds ~ ~ rHi(J,e + z/(J)r) 

as in |Tre02j . Nevertheless, the change we consider is better suited for points not 
extremely close to the double resonance. 

Recall that after switching to slow-fast variables. Hi does not depend on t. 
We have then the following changes 



e =0 - ^ log \y\ + X-\J)X'{J)xy + e £ Hi(J, 6 + u{J){s - X'^J) log | 2 /|))ds 

+ eTiy-\jy{J) (Hi(J, 6) - Hi(J, 6 - X-\J) log ||/|)) + 0*{e, xy) 
r^=J + eiy-\j) (Hi(J, 9) - Hi(J, 9 - p{J)X-\J) log bl)) + Oy,xy) 


T 


logbl 
■ A(J) 


+ 0*{e,xy) 


h =E{J) + X{J)xy + £Hi(J, 9) + Ol{e, xy) 


and 


le\ 


/e\ 

J 

= T 

V 

X 


T 

\y) 


\h) 


defined as 


(39) 


d y + + u{r])s)ds 

Jo 

+ X\y)X-\y)T {h - E{y) - 0) + eO*{e, h - E) 

J =V- £i^~\v) (Hi(r 7 , ^ + z/(r 7 )r) - Uy, Q) + eO*{e, h - E) 
X {h - E{t]) - eUy, 0) + h-E) 

y=ae^^^'^^ + 0*{e,h-E) 
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We call To the first order of this change of coordinates, dehned by 


/e\ 


(i\ 


/ ^ + u{p)t \ 

j 

= To 

p 


p 

X 

T 


^^e-^iv)r _ E{r^)) 

\y) 


\h) 


\ / 


(40) 


Theorem 3.4 can be rephrased as the following lemma. 


Lemma 8.1. Assume that the function wo in (20) satisfies 

^-i^i+a ^ \wo{ri*,h*,^*)\ < ce 

for some c > 0 and 1 > a > 0 independent of e. Then, the separatrix map SM.^ 
has the following form 


T] =1] 


+ ed^Q'^{ri,^,T) + eB^’‘^{ri,^,wo,T) 


r=e + /i" +3^ log 


K Wq 


X 


+ eB^’‘^{r],^,Wo,T,t) +0*{e) 


h* =h 

T* =T + t 


+ edrQ^ip, f, r) + r) + 


+ tlog 


A 


+ eB^’^{p,^,wo,T,t) +0*{e) 


a =a sgn wq 

where v, X, and k"" are functions of p, 

1 


Wq,t) = — {afip, ^ + z^r) - Hi{p, 0) 


- [iiiiv, ^ + /i'") - Hi(r/, ^ + - log 

= Hi(77,^ + z/r + /i'") -Ufip,^ + UT), 


K Wq 


X 


+ 


and 


B^’^{p,^,Wo,T,t) = fi{p,f,Wo,T)T + f2{p,^,Wo,T) log 
C, Wo, r, i) = gi{p, f, Wo, t)t + g2{p, C, Wo, r) log 


K Wq 


X 

kTWq 


X 


+ /3(h,C,wo,r,f) t 
+ fi'3(h,C,Wo,r,f) i 


for certain functions fi, pi which satisfy fi^pi = 0*{T). Therefore, the functions 
5^’°' satisfy 

B^’^,Bh’^ = 0*{1), B^’fiB^’^ = C>*(loge). 
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Recall that i/(0) = 0. Nevertheless, one can easily see that the function 
is well dehned even as z/ —)■ 0 since it has a well dehned limit. 

Proof. To prove this lemma, it is enough to compose the change of coordinates 


given in (39) with the map T as done in the proof of Lemma 7.2 One has also to 


take into account how to derive the splitting potentials 0^ from the dehnitions 
as explained in Lemma 7.3[ 


of the functions , 


□ 


A The separatrix map of the generalized Arnold 
example 


In this appendix we apply Theorems 3.1 and 3.4 to the generalized Arnold ex¬ 
ample ([^. The Arnold example presents several simplihcations and thus the 
formulas are considerably simpler than in the general case. For these models all 
the transformations and maps are C°°. The results presented in this appendix 
also apply to perturbations of the form 

:= P(J,p,exp(ig),exp(i(y9), exp(R)), 

where P is a real trigonometric polynomial of order N in ip and t, i.e. it has the 
form: 

4 - k2p + kst), 

\ki\<N,i=l,2,3 

for some N G 

The key point is that the generalized Arnold example has an unperturbed 
Hamiltonian Hq which does not present copuling between the “pendulum vari¬ 
ables” and the rotator. We hrst analyze the non-resonant regime as dehned in 


(16). Thus, we consider the resonant regime for any resonance k G 7V(Pi) U 
as dehned in dlT] ). 

For the generalized Arnold example, the frequency z/(p) dehned in ([^ satishes 
z/(p) = r]. The matrix A dehned in (10) is just 


A = 


0 1 
1 0 


therefore its positive eigenvalue A is A = 1 and is independent of I. The function 
g, dehned through the Moser normal form in Lemm a |4.1[ is independent of I and 
satishes g{r) = r + 02{r). Then, the function w in (|18[) is dehned by 


w = g 


-1 


h* 
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Analogously, since the Moser normal form is independent of J, so are the functions 
involved in the gluing map. Thus, the functions in ( |II| satisfy /i^ = 0. 
We have the following theorem. 

Theorem A.l. Fix /9 > 0 and 1 > a > 0. For e sujficiently small there exist 
c > 0 independent of e and a canonical system of coordinates {rj, h, r) such that 
in the non-resonant zone NR^ we have 

V = I + Ol (^e,Ho- 


P\ ( P 

^ = F + f ^ h = Hq Ol [ e, Ho —— 


where f denotes a function depending only on {I,p,q,e) and such that f = 0(e). 
In these coordinates the separatrix map has the following form. For any a G 
{ —,+} and {ri*,h*) such that 

^-i^i+a ^ \w{ri', h*)\ < ce, |r| < c~^, c < h*) \ P < c~^, 

the separatrix map {ri*,f*,h*,T*) = SAi(ri,ph,T) is defined implicitly as follows 


V* 

C 

h* 


7 ] — + 

^ + 

+ Ol(e log£) 
h- + 

T + t + 

+ Ol{e logs). 


+ Ogfyloge) 

dr,w{r]*, h*) [log \w{ri\ h*) \ + h*))] 

+ 0;(e) 

dhw{T]*, h*) [log \w{r]*, h*) \ + h*))] 


where M* and are C°° functions defined in Lemmas 6.1 and |7.^ respectively, 
and t is an integer satisfying &■ The functions M* are evaluated at h*, r) 

and satisfy 


+ o;{w), + o;{w), 

where G®" is the Melnikov potential defined by 


/ +00 

{Hi (T'^fy, pT + t),t-T) - Hi{r],^ + fy, 0, 0, f - r)) dt 

■OO 

(41) 

and r”" are the time parameterization of the pendulum separatrices, that is 

( 2a 

fy r) = {r],^ + r]T, 4 arctan(e'"^). 


cosh r 
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Proof. This theorem can be easily deduced from Theorem 3.1[ First, one needs 
to recall that the Moser normal form only depends on x and y (this is the reason 
for the better estimates for the function /). 

To deduce the formula of the Melnikov potentials 0°^ defined in (15) it is 
enough to recall that for the generalized Arnold model ([^, = 0. This implies 

that the contribution to (15) vanishes and therefore the splitting potential is 
just given by the Melnikov integral. Since Hq is uncoupled we also have = 0. 
Then, the Melnikov potential can be just written as one integral. □ 

Now we analyze the resonant regions. We consider the slow-fast variables 
(J, 6*) defined in Section]^ The function tco defined in ([^ just becomes 


Wo{v*,h*,e) = h* - 


iv* 




Theorem A.2. Fix /9>0, l>a>0 and k G Af{Hi) U For e 

sufficiently small there exist c > 0 independent of e and canonical coordinates 
[i], (^, h, r) such that in the resonant zone Res^ the following conditions hold: 

• the canonical form oj = dy A d^ + dh A dr; 


. r] = I + Ol{e,Ho-E{I)), i + yr = ip + f, h = Ho + Ol{e, Ho-E{I)), where 
f denotes a function depending only on {I,p,q,e) and such that f = 0(e). 

• In these coordinates SM.^ has the following form. For any a G , -|-} and 
{y*, h*) such that 




< \wo{y*,h*,C)\ < ce, 


r < c 


-1 


c < \wo{y*,h*,C)\e^ < c 


-1 


the separatrix map {y*, h*, t*) 
follows 

= SAi^{y,^,h,T) is defined 

implicitly as 

y* =y 

+ ed^e^{y,i,T) 

+ eB^’^{y,^,wo,T) 


c =e 

+ y\og \Hwo\ 

+ eB^’^{y,^,wo,T, t) 

+ 0*{e) 

h* =h 

+ edrQ'"{y,^,T) 



T* =T + t 

+ log Hwq 

+ eB^'^{y,^,T, wo,i) 

+ 0*{e) 

a* =a sgn wq 



where 0°" are 

the Melnikov potentials given m (41). If we 

define, z = 


{y,^,Wo,T) the functions 5'^’°' and are defined by 


i (Hi(r7,0 -Hi(r7,^ + 77log|K‘"'«;o|)) 
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and 


t) = fi{z)T + f 2 iz) log Ik'^woI + fsiz, t) t 
B^’^{z, i) = gi{z)T + g 2 iz) log |k'"wo| + i) t 


for certain C°° functions ft, gi satisfying fi,gi = 0*{1). 


Proof. To prove this theorem it is enough to use the results in Theorem 3.4 and 
take into account that the particular form of Hq implies that A = 1 and = 0. 
Then, one can easily see that B^’^^ and B^’^ have the form given in this theorem. 
Moreover, reasoning as in the proof of Theorem |A.l[ one can see that the splitting 
potential is given by formula 

□ 
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